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Abstract 

Significant advancement has been made in the last 
few years to identify noise sources in high speed jets 
via direct correlation measurements. In this technique 
turbulent fluctuations in the flow are correlated with 
far field acoustics signatures. In the 1970’s there was a 
surge of work using mostly intrusive probes, and a few 
using Laser Doppler Velocimetry, to measure turbulent 
fluctuations. The later experiments established “shear 
noise” as the primary source for the shallow angle 
noise. Various interpretations and criticisms from this 
time are described in the review. Recent progress in the 
molecular Rayleigh scattering based technique has 
provided a completely non-intrusive means of 
measuring density and velocity fluctuations. This has 
brought a renewed interest on correlation 

measurements. We have performed five different sets 
of experiments in single stream jets of different Mach 
number, temperature ratio and nozzle configurations. 
The present paper tides to summarize the correlation 
data from these works. 

I. Introduction 

The task of jet noise source identification is 
important for various industrial applications, most 
notably for aircraft noise abatement programs. In spite 
of the last half century of work, there is no unanimous 
agreement on the physics of jet noise emission. This 
has hampered development of physics-based prediction 
techniques. It is fair to say that a lot of discrepancies 
remain in predicting the low frequency, most energetic 
part of noise from a simple round nozzle using the 
existing models. The prediction techniques use 
simplified models of fairly complex two-point space- 
time statistics measured in real jets. The sources need 
to be multiplied by suitable Green’s function and 
integrated over the entire plume to recover the far field 
sound pressure fluctuations. 1 Experimental 
measurements of two-point statistics do not directly 


provide noise sources. Noise generation is 
fundamentally a coupling problem where energy from 
a small part of the turbulent fluctuations, that achieve 
sonic phase speed in the ambient, is ultimately radiated 
as sound. Separation of this tiny part from the large 
turbulent motion is indeed a difficult problem. It has 
become a common practice to use microphone arrays 
to determine sound sources solely from the far field 
measurements. Its advantage lies in its simplicity and 
the avoidance of measuring complex turbulent flow; 
the source being modeled simply as a distribution of 
monopoles. However, the noise sources are expected to 
have certain spatial coherence and preferred directivity 
which are yet to be correctly established. Besides, it is 
difficult to incorporate these parameters in the phased 
array processing. It can be argued that an independent 
means of source identification can neither be achieved 
by only turbulence measurement, nor by sole 
observations from the far field, but through a 
simultaneous measurement of the cause (turbulent 
fluctuations) and the effect (far field noise); this is the 
approach of the present study. 

By definition, correlation measurements require a 
priori knowledge of the far-field pressure fluctuations. 
Hence, it can not be used as inputs to the noise prediction 
codes. It is expected that there are three different usages 
of the correlation data. First, the validation of postulated 
source mechanisms in various jet noise theories. The 
experimental data may or may not follow the trends 
expected from certain theories. The second usage is 
simply as a diagnostic tool to identify changes in the 
noise sources caused by changing plume properties or by 
nozzle contouring. Third, the validation and debugging 
of computational aeroacoustics codes that strive to 
predict far field noise directly. 

Although first introduced more than 30 years ago by 
Siddon, 2 significant advancements in correlation 
measurements were made in the last five years with the 
advent of the molecular Rayleigh scattering technique. In 
the following results and observations from the earlier 
work are first presented. This is followed by an elaborate 
summary of results obtained over the last five years. 
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II. Prior Work 


The direct correlation between the cause and effect 
was originally proposed by Siddon 2 (also see Rackl 3 ) to 
locate sound sources. A surge of enthusiasm appeared 
during the 1970’ s when many experimentalists took up 
the method. The primary objective of these experiments 
was to validate various theories such as the self and shear 
noise theory of Ribner, 4 fluid dilation model of 
Meecham & Ford. 5 and coherent pressure fluctuations 
model of Michalke & Fuchs. 6 The first one obtained 
most of the attention, where far field sound pressure 
fluctuations were correlated to the turbulent velocity 
fluctuations inside the jet. 7-10 The last two required 
correlations with the pressure fluctuations 1 1-12 inside the 
jet. The noise generation theories are still discussed, but 
significant doubts have been expressed in many of the 
experimental results. Usually a microphone with a nose 
cone is inserted into the flow to measure pressure 
fluctuations. Very high correlation values ~0.4 were 
measured between the flow and far field pressure 
fluctuations; 12 however, it is difficult to separate out the 
true correlation values from that caused by probe- 
interference. The same is true for velocity measurements 
by hot-wire probes. 7 ’ 8 Richarz 10 writes, ’’Although the 
noise radiation from a hot-wire probe is completely 
masked by jet noise, it may contribute to the measured 
jet flow-jet noise correlation nevertheless. A little 
thought will show that the unsteady forces acting upon 
the probe are proportional to the velocities measured by 
it. Therefore, the probe noise and velocity signal are well 
correlated.” With this realization a small number of the 
early researchers, most notably Schaffar 9 and Richarz, 10 
resorted to non-in vasive Laser Doppler Velocimetry 
(LDV) for flow-fluctuation measurements. The rest of 
this section will discuss some of their results. 

The theoretical basis for the correlation measurements 
are from Proudman’s version of the Lighthill’s acoustic 
analogy equation. The far field pressure fluctuations p'at 
the observer position X f is due to the turbulent 
momentum flux pu r u r along the direction of the observer 
(p is the gas density, and u r is the velocity vector along 
the observer direction, figure 1): 


P(X f ;t) = 


4ita5r 


J^[pu r u r (X s ;t')]d 3 X s (1) 

\7 & l 


Where X s is location inside the jet, r = X r X s is radial 
distance, a 0 is the ambient sound speed, and the 


integration is carried out at a retarded time 
t 1 = t - To, T 0 = |Xf-X s |/a 0 The above description neglects 
viscous effects and the contribution from the “second 
source,” i.e., temperature fluctuations. To obtain 
acoustic intensity an auto-correlation function needs to 
be worked out. In the “causality” formulation this is 
accomplished by multiplying the source integral with 
the far field sound pressure fluctuations: 


p ,2 )(X f ;x) = 


f^fpu r u r (X s ;t'+T)]d 3 X s [p'(X f ;t)]' 

In afj r y 3t 

= 777^I I '-"' (X ‘’ X ' ;,| ' |I>X ' (2) 

4n r a 0 ydx 

R P u r ii r ;p'( X s>X f ;x) = (pu r u r (X s ;t + x-x 0 )p'(X f ;t)) 


The < > bracketed expressions imply time averages. 
Such time averages are also expressed with an 
additional semi-colon, e.g. <pu r u r ; p'>, in the text. The 
above equation says that the auto-correlation of the 
sound pressure fluctuations can be obtained from the 
cross correlation of the fluctuating source terms and 
the radiated noise. 

In the past, density fluctuations have been neglected 
from the source term; although it introduces error, 
especially for heated and higher Mach number jets. (In 
the primary plume of a typical turbofan engine 
p 1 / p ~ 0.4). Ribner also decomposed instantaneous 
velocity into time-averaged and fluctuating parts: 

u r =u r +Up. This splits the flow-sound cross- 
correlation into two parts: 

R pu r u r ;p , ( x s 5 X f;x) = 2pu r (u r (X s ;t + x-xo) p (Xf ; t)^ 
+ P (u' r u'r (X s ;t + x-x 0 )p'(X f ;t)) (3) 

The first term on the left is the ‘shear noise’ source and 
the second is the ‘self noise’ source. The relative 
importance of these terms in noise radiation at different 
polar angles is still under debate. One important 
observation of Schaffar 9 and Seiner and Reethof 8 is 
that in the shallow angle to the jet axis the ‘shear 
noise’ exceeds ‘self noise’ by at least 17dB. 

In both LDV and hot-wire based measurements only 
the axial component of velocity u was measured; the 
component along the observer direction was 
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determined by multiplying with suitable direction 
cosines. The low speed jet, Mj = 0.3, used by Richarz 10 
and Seiner & Reethof 8 produced extremely small 
correlation coefficients. The latter had to use 
narrowband processing to get any signal at all. Typical 
values of the normalized correlation <u 7 ; p / >/u rins p rms 
was 0.005 (subscript rms stands for root-mean-square) 
in Richarz’s experiment. Schaffar, on the other hand, 
used a higher speed Mj = 0.98 jet, where maximum 
correlation value increased to 0.05. This dependency of 
the correlation coefficients on jet velocity was 
confirmed later from the Rayleigh scattering based 
measurements. 13 ' 14 The ultimate success of the prior 
causality work appeared when Schaffar calculated the 
volume integral in equation (2) from a plume survey. 
For correlations measured from a fixed microphone, 
positioned at shallow 160° and 150° angles, time 
derivative and integration following equation (2) 
reproduced the microphone auto-spectrum (fig. 2). 
Although a large number of approximations had to be 
made, yet the result showed that axial velocity 
fluctuations u 7 (not uV the self-noise source) are 
responsible for noise at shallow angles. Schaffar failed 
to measure any correlation with either u 7 or u 72 terms 
for microphone polar angle 0 P > 135°. Such a steep 
drop in correlation coefficients is also confirmed in 
recent measurements. 

Publications on correlation study for jet noise 
became sparse at the beginning of 1980s. Perhaps, the 
realization of the intrusiveness of the easily available 
probes, criticisms on the interpretation of the data, and 
receding interest on the jet noise problem contributed 
to the decline. 

III. Recent Work 

The overarching interest of the prior work was to 
validate different noise generation theories via 
correlation measurements. Efforts to validate the 
causality principle of equation (2) using recent data 
produced confusing results. At certain microphone polar 
angles no correlation can be measured; therefore the 
causality relationship cannot be applied. Yet at certain 
other polar angles correlations were so high that a unit 
volume at the end of the potential core is found to 
produce more noise than the entire jet. 13 Therefore, the 
primary focus of the recent work was to find the nature 
of the correlation coefficients from different plume 
conditions and nozzle configurations without definite 
association with any theory. 


The correlation technique has been revisited in the 
last few years with the aid of a new molecular 
Rayleigh scattering based instrumentation that is 
capable of measuring unsteady density and velocity 
fluctuations in heated, high speed flows. 14 ' 15 
Fluctuations occurring over a frequency range of 0 to 
50 kHz have been measured. The technique is based on 
laser light scattering from the gas molecules present in 
air. Since neither any seed particles nor any intrusive 
probes are used, the technique is truly non-invasive. 
Implementation of the technique requires special 
attention to the cleanliness of the air stream, isolation of 
the sensitive optical components from jet noise, usage of 
Fabry-Perot interferometer, low-level light measuring 
electronics and minimization of the impact of electronic 
shot noise during signal processing. 

Molecular Rayleigh scattering technique. — To limit 
the length of this paper detail description of the 
instrumentation is avoided; that can be found in 
references 13 to 17. An advantage of the Rayleigh 
scattering technique is that in addition to the 
fluctuations in any velocity components, fluctuations 
in scalar variable, such as density (and in theory 
temperature), are measurable. So far, advancement is 
made in simultaneous measurement of density and one 
component of velocity. The central element is a 
narrow, continuous wave laser beam which is passed 
through the jet plume that is cleansed of dust particles. 
Scattered light from a small length of the beam was 
collected using a set of optics and then analyzed to find 
the flow properties. The probe volume diameter was 
about 0.15 mm, and length varied from 0.5 to 1 mm 
between setups. For a fixed mixture of gases, such as 
air, the power of the scattered light is directly 
dependent on local density. 14 To measure the time 
variation of density the variation of the scattered light 
intensity needed to be monitored. This was 
accomplished by using a photo-multiplier tube and 
photo-electron counting electronics. To measure 
velocity, the Doppler shift of the scattered light from 
the incident light needed to be measured. An optical 
spectrum analysis of the scattered light was performed 
by passing the scattered light through a Fabry-Perot 
interferometer. 15 The Fabry-Perot was operated in an 
imaging mode and the fringe formed at the image 
plane was monitored. Only a fraction of the free- 
spectral range was used, and all frequency shifts were 
measured from a single fringe. At first, the 
interferometer was placed into a control loop that 
stabilized the position of the reference fringe created 
by the unshifted incident light. When Doppler shifted 
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light was passed through the interferometer the fringe 
was formed at a different radius. (It’s width is also 
increased due to thermal broadening). To measure the 
time variation of velocity, the time variation of the 
shifted fringe needed to be monitored. This was 
achieved by a pair of concentric mirrors that split the 
image formed by the interferometer into two parts. The 
light intensities from the two parts of the fringe were 
measured by a pair of PMT. An in-situ calibration 
process was employed where the ratio of the photo- 
electron counts from the two PMT were measured for 
different flow velocities. In unknown flows, photo- 
electron counting was performed in contiguous time 
bins, which upon multiplication by the calibration 
constants provided a time history of velocity 
variations. The scattering geometry was setup to 
measure either the axial component of velocity u or the 
radial component v. 

All data presented in this paper are from two 
different single-stream jet facilities at the NASA Glenn 
Research Center. Table I summarizes the operating 
conditions and nozzle configurations. Figure 3(a) 
shows a photograph of an unheated jet facility where 
one convergent nozzle, two C-D nozzles of design 
Mach number Mj = 1.4 & 1.8, one 16-lobed 
rectangular nozzle and one 4-tabbed circular nozzle 
were tested. 13 ' 14 ' 16 Table I describes jet velocity in 
terms of the fluid dynamic Mach number Mj as well as 
the acoustic Mach number M a . Since noise radiation is 
a coupling process where flow fluctuations acquire 
acoustic speed in the ambient, acoustic Mach number 
is believed to be of more relevance. 

In addition to single microphones, a microphone 
phased array was also employed; 18 Figure 3(a) shows 
the setup with a 32 channel array of electret 
microphones mounted on a semi-cylindrical “bucket” 
located just below the laser probe volume. The array 
was concentric with the jet cross-section. It was also 
traversed with the laser probe volume which was 
located at the center of the array. The semi-circular 
array had a radius of 9.25” and the microphones were 
arranged in two concentric ellipses. Figure 3(b) shows 
the Rayleigh setup around the heated jet facility. In 
both facilities the primary air stream was filtered to 
remove dust particles. In addition a low speed 
(~15 m/s), clean, co-flowing stream was created 
around the primary jet to block entraining of the dust- 
laden ambient air. A hydrogen combustor was used in 
the heated facility to avoid soot particles. The bulk of 
the correlation data are from single microphones 
mounted on a 50D (D: nozzle exit diameter) arc for the 


unheated jet and 100D for the heated jet facility. The 
laser probe volume was traversed on an x-r plane 
(x: axial, r: radial directions) containing the jet 
centerline. 

Comparison of correlations from different flow 
variables. — The unheated M a = 0.91, 1.18 and 1.4 jets 
were investigated using two different Rayleigh setup 
which measured either p and u, or p and v 
simultaneously (u is the axial and v is the radial 
component of velocity). Figure 4 presents normalized 
<p; p ; >, <puu; p> and <pvv; p > cross-correlations from 
150° and 90° microphone polar angles. The puu 
product was created by multiplying density p(t) and 
velocity u(t) time -traces. The cross-correlation values 
were calculated via Fourier transform, that is, the 
cross-spectrum P puu; p/ and P p: p/ were calculated first, 
and then inverse transforms were taken to return to the 
time domain. A significant problem with the Rayleigh 
time series data arose from the electronic shot noise 
which added a large random white noise to the velocity 
and density signals. The cross-correlation process 
minimized contamination from shot-noise. However, 
the cross-correlation data needed to be normalized by 
standard deviations of puu and p fluctuations which 
were also contaminated by the shot noise. To reduce 
this contribution, the standard deviations were 
calculated from individual spectra from which the shot 
noise floor was subtracted. The latter was estimated 
from the highest frequency bin where energy from the 
turbulent fluctuations was deemed to be minimal. It is 
estimated that the noise removal process created 
±10 percent error in the correlation data. 

The cross-correlation data of figure 4 show sharp rise 
at a time delay which slightly exceeds the time needed 
for sound waves to travel from the laser probe location 
to that of the microphone with ambient speed, 
To = (X f - X s )/a 0 . The laser probe was located at the 
centerline and farther downstream from the end of the 
potential core: a region found to provide the largest 
correlation. There are multiple interesting observations 
that can be made from figure 4. First, air density 
fluctuations show as good a correlation as puu; even 
the time variations of the correlation data are similar. 
The same was found to be the case for the correlation 
with the axial velocity fluctuations <u; p >. Therefore, it 
can be said that the noise sources can be expressed in 
many different variables. There are, however; other 
variables which show far poorer correlation. This is the 
second observation. Unlike puu, pvv show significantly 
lower correlation with the far field noise. That the 
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<pvv; p'> correlation with 90° microphone is 
immeasurably small is of particular interest. The 
causality relation of equation (2) dictates that 

( pvv ; \) '"j correlation as the source for the 90° noise. 

Therefore it had been expected that <pvv; p'> 
correlation from 90° microphone would be significant, 
which was found to be incorrect. Similarly causality 
relation shows that pvv fluctuations has little to do 
with 150° noise radiation, yet microphone signal from 
this angle show significant ( pvv ; p ''j correlation above 

the noise floor. Third, sound pressure fluctuations at 
90° to the jet axis correlate poorly with any flow 
variables. An exception was found in <p; 
p>corrclation in highly heated jets which will be 
discussed latter. Fourth, the time duration Ax over 
which correlation changes from zero to negative to 
positive and back to zero is significantly long. For the 
particular M a = 0.95 jet, Ax is estimated as 1.2 msec. 
This provides a measure of coherence time of the 
dynamic process responsible for sound radiation. 
Invoking Taylor’s hypothesis, the corresponding 
coherence length scale is l co her = Ax U c , where U c is the 
convective speed. Assuming U c =0.7Uj the coherence 
length koher ~ 10D. In other words, the longest source 
contributing towards correlation is, on an average, 10 
jet diameters long. Therefore, the noise source 
responsible for sound radiation, at the shallow 150° 
angle, is indeed “non-compact.” 

Since many of the earlier studies revolved around 
self and shear noise sources, attempts were made in 
reference 13 to revisit this division, albeit in a more 
general way. Unlike the earlier efforts density 
fluctuations are now measurable; hence there appear 
additional shear noise term from p 1 fluctuation, and self 
noise terms from pV and pV if fluctuations: 

p = p + p', u = u + u' 

(puu;p') = (2puu';p') + (p'uu;p') + (4) 

(2 up'u p') + (p u 'u p') + (p'u'u p') 

The five fluctuating terms in puu decomposition were 
constructed from the measured p(t), u(t) time traces 
and individually correlated with the sound pressure 
fluctuations. Table II shows relative importance of 
various terms calculated by uniform normalization. 
This table confirms that the terms containing first order 
fluctuation, attributed as shear noise source, provides 


far higher correlation than the second and third order 
fluctuations, the self noise sources. 

Acoustic Mach number M a and polar angle 0 P 
dependence". — The largest value of either <puu; p'> 
or <p; p'> correlation is found to be a strong function 
of the acoustic Mach number. Figure 5 presents plots 
of the maximum <p; p'> correlation measured with the 
sound pressure fluctuations from microphones kept at 
different polar angles. Similar dependency in the <puu; 
p'> correlations is presented in reference 13. The 
normalized correlation coefficients can be interpreted 
as the efficiency of turbulent fluctuations present at the 
laser probe volume to radiate to the microphone 
location. Clearly, such radiation efficiency improves 
with the Acoustic Mach number. The most significant 
improvement is seen for fluctuations along the lip line. 
They are mostly non-radiative in the subsonic 
M a = 0.6, yet become efficient radiators in the 
supersonic M a = 1.48 case. Polar-angle-wise, sound 
pressure fluctuations from the shallowest angle of 150° 
show the highest correlation. Correlation values from 
the 140° and 130° microphones are also significant. 
There occur sharp drops around 120° and especially 
for the lower M a jets. For the higher M a jets some 
correlations are measured even for the 90° microphone 
when the laser probe is placed around the end of the 
potential core. It is believed that the present correlation 
study, for the most part, is identifying turbulent 
fluctuations that radiate close to the jet axis. 

Effect of heating . — To explore the effect of heating 
on noise sources, experiments were conducted in a 2” 
diameter heated jet facility mentioned earlier. (The 
unheated jet data in figure 6(e) is from the 1” facility). 
To isolate the effect of heating the plume velocity (and 
therefore the acoustic Mach number M a ) was kept 
constant, while the jet static temperature was increased 
by hydrogen burning combustor chamber. Figure 6 
presents data from a large range of temperature ratios 
and demonstrates a direct relationship between the plume 
temperature and correlation coefficients: an increase in 
one is reflected similarly in the other. This increase, once 
again, is most reflected in the peripheral shear layer 
which hardly radiates in the unheated subsonic jet, yet 
becomes increasingly efficient as T r is increased. This 


* In references 13, 14, and 16 and prior work in the 1970s 
microphone polar angles were measured from the jet flow 
direction, while for the present document polar angles are 
presented from the flight direction: 150° in the present paper 
is equivalent to 30° in the earlier. 
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increase in radiation efficiency may be due to the second 
source associated with temperature fluctuations in 
various acoustic analogy equations. 

In order to demonstrate the effect of heating, 
reference 17 presented the coherence spectra of the 
correlation values. Unlike the time domain correlation 
values presented in figure 6, coherence spectra are 
heavily affected by electronic shot noise. This was 
particularly true for the T r = 1.0 jet where the plume and 
the ambient had nearly the same time averaged density. 
These led to very small level of fluctuating density and 
ultimately nearly zero coherence. Some of the 
conclusion drawn in that paper based on the coherence 
spectra may not be correct. 

Frequency dependence and distributions of noise 

sources. — So far data from time-domain correlation 
study were discussed. Such data do not show the 
frequency dependence of correlation. Towards this 


goal, coherence spectrum F^/ /(fj) was calculated 


from the time series of the photoelectron counts N 
(which multiplied by calibration constants provide air 
density) and microphone pressure fluctuations: 


r 2 f P NV (f l } 

rp/p/(fl) = P N^N^ f l) P pV (f l) (5) 

where, P N , p , (fj) = (f n , (fj ) . F* (f , ) ) 


Here, F / represents Fourier transform of time- 

series of photo-electron counts, F*/ is the complex 

conjugate of the Fourier transform of the microphone 
signal, ? N / N / and P / / arc corresponding power 

spectra, and n is the total number of samples used for 
the calculation. The center frequencies fi were then 
converted into Strouhal number (St = fD/Uj, D: jet 
diameter, Uj = jet exit velocity). One important 
disadvantage of the coherence spectra is that the effect 
of shot-noise can not be entirely removed. The cross- 
spectrum mostly rejects this noise, but the power 
spectrum of photo-electron count contains contribution 
from both the desired density fluctuations and the 
undesirable shot noise: 



calibration cons tan t 


+ a 


2 

shot noise 


( 6 ) 


Where, cr p and a 2 s h 0 t noise represents mean square of 
respectively, density and shot noise fluctuations. In the 
unheated low speed jets the later can be much higher 
than the former, resulting in a large bias towards the 
lower value in coherence spectrum. The situation 
worsens for a temperature balanced jet, discussed 
earlier, where coherence spectrum falls mostly in the 
noise floor. Nevertheless, the coherence spectrum 
r 2 pp (St) provides a measure of linear dependency 
between turbulent density fluctuations and the far field 
sound pressure fluctuations. In the absence of shot noise, 
a coherence value of unity implies a perfectly linear 
cause and effect relation while a value of zero implies no 
correlation. Convergence errors and electronic shot noise 
never allowed attainment of the zero value even in 
quiescent flows; instead a noise floor around r 2 pp - - 
0.002 was reached. When the coherence level is above 
the noise floor, it can be said that some part of density 
fluctuations from the laser probe location is creating 
sound pressure fluctuations at the microphone location; 
in other words, the probe location is a sound source. 

The spatial distribution of coherence values at different 
Strouhal frequencies for two different Mach number jets 
is shown in figure 7. For this figure a large number of 
correlation data were measured by moving the laser 
probe from point to point in the flow, while keeping the 
microphone location fixed. Later on, coherence values at 
a desired Strouhal frequency were isolated for all 
measurement stations and plotted using indicated color 
scale. Note that the color scale is such that the minimum 
value (green for T = 0.005) is above the noise floor 
(r = 0.002). Therefore regions of no or very low 
correlation do not show any color. Since no correlation 
was measured in Mach 0.95 jet for St > 0.4, coherence 
distributions at St = .5, 1 and 1.5 are not plotted. A clear 
distinction in measurable sound sources in subsonic and 
supersonic jets is visible in this figure. Density 
fluctuations causing sound generation up to St = 1.5 
could be determined only in the supersonic case; the 
higher the frequency the closer the source to the nozzle 
exit. Additionally, the peripheral shear layer around the 
potential core is found to produce no correlation in the 
subsonic jet, while significant correlation at high 
frequency is measured in the Mach wave emitting 
supersonic jet. The lower frequency source in both jets is 
the strongest along centerline, after the end of potential 
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core. Such sources weaken rapidly in the radial direction 
and very slowly in the downstream direction. 

Effect of nozzle configuration. — A utility of the 
correlation method is that it can be used to identify 
changes in the source distribution made by nozzle 
contouring or other noise abatement schemes. A 
comparison of coherence spectra from different nozzle 
configurations demonstrate this utility (fig. 8). A 
baseline axisymmetric, a 16 lobe rectangular, and a 
4-tab axisymmetric nozzles were used to produce 
unheated Mj = 0.95 jets. 16 All three nozzles had 
equivalent diameter D eqv of 1”. The far field noise 
spectra presented in part (d) of this figure shows that 
within the Strouhal frequency range of 0 < St < 0.6 the 
axisymmetric configuration is the loudest and the 
lobed one is the quietest. The scenario reverses at the 
higher frequency end. The correlation data in the low 
frequency side reflects the same trend as seen in the 
noise spectra. The air density fluctuations present in 
the axisymmetric configuration is the most effective in 
sound radiation while the lobed nozzle is the least; the 
tabbed nozzle falls in between. Note that the coherence 
data from the tab-nozzle and the lobe-nozzle are shown 
for axial stations closer to the nozzle exit; this is to 
reflect the faster mixing and shorter potential core 
length which brings noise sources closer to the exit. 
Detailed flow data from these jets were presented in 
reference 16. Since, no correlation can be measured for 
St > 0.6 nothing can be said about the high frequency 
noise sources. 

An interesting double hump behavior in the 
coherence data from some stations (x/D eqv = 3, 4 
figure 8(c)) in the tabbed jet plume led to further radial 
surveys shown in figure 9. The coherence data from all 
radial stations at x/D = 3 show the same double hump 
behavior indicating the existence of two distinct noise 
sources. It is conjectured that the two sources may be 
due to the splitting of the jet into a bigger core jet and 
4 smaller jetlets outlined in the figure 9(b). The lower 
frequency part may be attributed to the global length 
scale from the jet diameter and the higher frequency to 
the interactions of the small jetlets. Notably, Tam & 
Zaman 19 has proposed a similar jetlet model to explain 
increased high frequency noise from tabbed nozzle. 

Microphone phased array. — Currently the most 
prevalent experimental tool to identify noise source is a 
microphone phased array. The distribution of the noise 
sources obtained by the present correlation technique 
and by phased array 20 carries some similarity and some 
dissimilarity. For example, both determine the end of 
the potential core as the strongest source for the low 


frequency noise. Phased array can detect higher 
frequency noise sources in any jet; the correlation 
method mostly fails to do so except for the high Mach 
number conditions. For these instances the progressive 
upstream shift of the noise source with increase in 
frequency is a similar result from both techniques. The 
important differences appear in discriminating noise 
sources radiating to different polar angles, and spatial 
resolution in source identification. The correlation 
technique is deemed superior in these regards. Usually 
phased arrays cover a large polar angle and the 
inherent monopole directivity assumption forbids 
discrimination based on the noise directivity. The large 
spatial extent of the instrument function translates into 
poor spatial resolutions in the phased array 
measurements. For example, the different behavior of 
the turbulent fluctuations along the lip line and the 
centerline found in the correlation measurements are 
difficult to reproduce using phased arrays. There are 
ongoing efforts to improve spatial resolution via 
deconvoluting the instrument function from the 
measured distribution 18 and by employing caged 
microphone arrays; still resolution obtainable by the 
point measurement process in the correlation technique 
is superior. Nevertheless, phased array technique 
avoids turbulence measurement which explains its 
omnipresence. 

It was conjectured that the lack of correlation in the 
high frequency range may lie with the usage of a single 
far field microphone. Since noise sources are 
distributed over the entire plume, the contribution from 
the small probe volume was overwhelmed by radiation 
from large number of distributed sources. A 
microphone phased array, on the other hand, is capable 
of extracting sound pressure fluctuations from a 
localized region of the jet. Therefore, a combined 
phased-array and Rayleigh scattering experiment was 
undertaken. 18 Figure 3(a) shows a photograph of the 
semi-cylindrical array of 32 microphones mounted on 
the jet facility. Note that the array was designed to 
cover a large azimuthal angle and a small polar angle, 
which made a narrower instrument function in the 
radial direction and a larger one in the axial direction. 
In other words, for the same Strouhal frequency the 
array was able to localize sources more precisely in the 
radial direction than in the axial direction. 
Experimental data from the combined phased array and 
microphone measurements can be found in 
reference 18. The large number of microphones used in 
the test allowed examination of the azimuthal variation 
of the correlation coefficients and the nature of the 
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acoustic phase fronts emitted by the jet. Only these 
parts are described in the following. 

Figures 10(a) and (b) allow a closer examination of 
the correlation data from different microphones located 
inside the phased array. Note that the microphones 
used in the part (a): number 0, 23, 29, and 10, lie 
parallel to the jet axis with 0 being the most 
downstream and 10 as the most upstream. The 
separation between these two is 6. 2D. The laser probe 
volume lies at the center of the pattern. Microphone 
used in figure 10(b): number 5, 20, 23, 26, and 15, lie 
at fixed axial position (x/D = 8) but at different 
azimuthal angles (microphone 23 is slightly upstream). 
Figure 10(b) is easier to explain: since the correlation 
data are all in phase, the sound source associated with 
the density fluctuations at the probe location is 
axisymmetric in nature. The different shapes of the 
correlation profile in the figure 10(a) perhaps indicate 
axial distortions of the phase front radiated from the 
probe location. 

Figure 10(c) shows correlation with the beam- 

formed signal p b f: - — — . The beam-formed 

Prms (Pbf ) rms 

time series was calculated by applying estimated 
propagation delay T; from the probe volume to the 
individual microphones: 

1 31 

P bf (0=— XPi( t+T i) (7) 

i=0 

The correct estimate of the time delay lies at the heart 
of the beam-steering process employed in the phased 
array processing. The non-compactness of the jet noise 
source, as evident from the long time duration of the 
correlation data makes it difficult to choose a single 
time delay. By increasing the time delay one may 
intend to steer the array to a further location, but that is 
also equivalent to looking at a different phase of the 
acoustic front arriving from a nearby point. This 
exemplifies the central problem of beam forming when 
the noise source is non-compact. For the effort of 
improving correlation via beam-forming, it was 
decided to use the simplest free-space propagation time 
for time delay estimates. Figure 10(c) shows an 
improvement of correlations by a factor of two. Data 
obtained from some other locations, for example the 
peripheral shear layer close to the nozzle exit, 
however, did not show much improvement. A 
frequency analysis, showed the same absence of 


correlation for St>0.6 found in the single microphone 
data. 

IV. Interpretation of Correlation Data 

The large value of correlation measured from certain 
locations in the jet and microphone polar angles and its 
near absence from the other are indicative of 
interesting noise source behavior. Attempts to 
associate correlation data with various noise generation 
theories are bound to generate debates. The following 
describes some criticisms and plausible explanations. 

Criticism of the causality approach’s. — The primary 
criticism of the causality approach is that the existence of 
correlation is not sufficient proof that the measured 
turbulent fluctuations are the sources of far field noise. 
Both may be the result of an unidentified forcing 
function. It is interesting to read criticism of Professor 
J.E. Ffowcs Williams and rebuttal from Professor T.E. 
Siddon in a discussion session, at the end of the 1973 
AGARD meeting. 21 Ffowcs Williams: “...if you were to 
measure the sound produced by my voice at the lips of 
Prof. Kiichemann and at the ears of somebody beyond 
him you would find that the correlation between the 
pressures was perfect. Yet is Prof. Kiichemann talking or 
am I?” In rebuttal Prof. Siddon replied,” . . .correlations, 
if normalized by the individual root mean square values 
of the partner variables, will give an erroneous 
impression of the source distribution. For example, two 
microphones both in the far field of, but on a radial line 
from a complex source will give a maximum correlation 
coefficient of unity, when normalized in the foregoing 
manner. In contradistinction the un-normalized causality 
functions... yield legitimate distribution functions which 
quickly drop to zero outside the source region..” Perhaps, 
“quickly drop to zero” should be replaced by: quickly 
drop to very small value. Nevertheless, the value of 
Prof. Siddon’ s rebuttal becomes apparent by examining 
the un-normalized values of <puu; \)> correlations. The 
largest fluctuation present in a jet is due to the fluctuating 
axial momentum flux puu, and this term is seen to be 
well correlated with the shallow angle noise. Radial 
variation of un-normalized correlation, measured by a 
moving probe and a fixed microphone, is show in 
figure 11. This figure shows that the un-normalized 
<puu; p 7 > drops by a factor of 100 when the laser probe 
is moved from the jet centerline to the edge of the 
turbulent plume. There exist p, u and p fluctuations in 
the free-field outside the jet boundary. But the absolute 
values of these fluctuations are so small that the laser 
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probe is unable to detect their presence. The difficulty of 
separating the source fluctuations from that of a passing 
sound wave is that the fluctuations in the same variables 
appear in both. Therefore, the physical expectation of the 
strongest fluctuations as the source and the weaker 
fluctuations as the effect ultimately differentiates the 
two. 

Turbulence as a sum of random eddies. — 

Correlation coefficients can be interpreted as the 
radiation efficiency of the fluctuations present at the 
probe volume. In the past efforts were made to 
associate the fluctuations with individual turbulent 
eddies. 7 ' 8 This stem from a picture of turbulence which 
is well described by Lee & Ribner: 7 “A turbulent jet 
may be regarded as an assemblage of ‘eddies’ 
(correlation volumes) which radiate sound 
independently and incoherently; the pressure 
contribution from each adds in the mean square at 
points in the far field.” Based on their measured 
correlation of 0.002 Lee and Ribner estimated that 
about 2500 eddies of equal strength radiates to the far 
field. This was used to determine sound emission per 
unit length of the jet. An examination of the earlier 
data shows that the measurable correlation was found 
only in the low Strouhal frequency range 0 to 0.5. It is 
now known that, over this frequency range, turbulence 
as incoherent eddies provides an incorrect description 
of reality. Therefore, the idea of calculating number of 
eddies generating sound is no longer persuasive. 

Large-scale versus fine-scale noise source. — 
Various aspects of the correlation data are better 
explained from a structural description of jet 
turbulence (Tam et al., 22 Michalke 23 ). The sound 
pressure fluctuations measured by a far field 
microphone are due to a sum of the radiation from all 
turbulent eddies distributed in the entire jet plume. 
However, based on the extent of spatial coherence, the 
eddies can be broadly divided into large structures with 
coherence length scales of the order of a jet diameter, 
and fine structures of much smaller spatial coherence. 
Figure 12 provides a schematic description of the 
situation. It is believed that the present experiment 
picks up contribution from the large structures and is 
unable to determine contribution from fine scale, small 
eddies. Although the current technique measures 
correlation from a single point in the plume, the 
correlation values are extremely high for shallow 
angles. The radiation from various spatial locations of 
large structures is expected to be phase related; 
therefore, a single point correlation reflects 
contribution from the entire correlated source. 


Theoretical calculations of Michalke 23 show that the 
higher the spatial correlation of turbulence, the 
narrower the radiation angle. Additionally, frequency 
analysis of cross-correlation data demonstrates that the 
measured correlations are from low Strouhal frequency 
fluctuations. These are characteristics of large 
organized structures which radiate primarily in the 
downstream direction close to the jet axis. The quick 
drop in correlation with an increase of microphone 
polar angle, seen in figures 5 and 6, is reflective of the 
inefficiency of the large scale structures to radiate at 
90° or upstream angles. According to the model of 
Tam et al. 22 fine-scale structures with little spatial 
coherence are primarily responsible for sound radiation 
at 90° and lower polar angles. Such eddies radiate 
more omni-directionally; the net sound field at a far 
field point is a sum of the contributions from many 
such eddies with random phase relationships. Hence, 
correlation from a single measurement point in the 
plume is expected to be very small, below the 
experimental noise floor. This perhaps explains 
inability of the present point measurement scheme to 
identify noise sources for the 90° radiation. 

V. Concluding Comments 

The correlation method to identify noise sources has 
remained an exclusively experimental effort. 
Computational aero-acoustics did not exist in the 1970’s. 
It will be of interest to see modern computational efforts 
to replicate some of the observations made from the 
experiment. The experimental data can also be used to 
determine the resolvable frequency range in the large 
eddy simulation and the usefulness of direct numerical 
simulation for predicting noise from large Reynolds 
number jets. 

Experimentally the lack of correlation for St > 0.5 
(St > 1 for Ma > 1 jets) has remained a challenge. 
Further innovative experiments are required for the 
future progress. The computational approaches can also 
be used to diagnose the underlying reasons. 
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Table I**. — Operating condition, nozzle geometry etc. 


Reference 

Acoustic 
Mach no, 
M a =U/a 0 

Jet Mach no., 
Mj = Uj/a, 

Static temp ratio, 

T, 

Nozzle config. 

Comments 

13, 14 

0.911 

0.95 

0.853 

Round 

Also investigated using 
mic phased array (ref. 1 8) 

13, 14 

1.18 

1.4 

0.725 

Round 

C-D nozzle 

13. 14 

1.4 

1.8 

0.614 

Round 

C-D nozzle 

13, 14 

0.58, 0.67, 
0.77, 0.84 

0.6, 0.7, 0.8, 0.9 

0.93,0.91, 0.886, 
0.86 

Round 

Limited data 

17 

0.9 

0.982, 0.9, 0.754, 
0.671, 0.604 

0.84, 1.0, 1.429, 
1.818, 2.273,2.7 

Round 

Cases with fixed M a but different 
TR 

17 

0.6, 0.9, 1.48 

0.402, 0.604, 0.993 

2.273 

Round 

Cases with fixed TR but 
different M a 

16 


0.95 

0.85 

Round, tabbed 

Asymmetric nozzle 

16 


0.95 

0.85 

Rectangular 
lobe nozzle 

Rectangular nozzle 


**Typical operating condition of a modern turbofan engine: Primary flow: M a = 1.5, Tr = 2.27; Fan flow: 
M a = 0.98, Tr= 1.01. 


Table 11 . — Peak correlation coefficients between far field sound pressure fluctuations and various 
Reynolds decomposed terms of puu: 


Mach 

no. 

Probe 

x/D 

Probe 

r/D 

Mic 

ang 

(puu;p^ 

^2p UU 7 ; p 7 ^ 

(p^mp^ 

(^2up / u / ;p / ^ 

(puu)rmsPms 

(puu)rmsPms 

(p uu )rms Prms 

(p uu )rms Prms 

1.8 

12 

0 

150° 

0.2207 

0.1635 

0.0636 

-0.00214 

1.8 

12 

0 

SO 

O 

O 

0.0056 

0.0033 

0.0023 

Noise 

1.4 

10 

0 

150° 

0.1919 

0.1344 

0.0624 

-0.0021 

0.95 

10 

0 

150° 

0.0682 

0.0592 

0.01 

-0.0002 


^puV;p^ 

^p / u / u / ;p / ^ 

(puu ) rm s Prms 

(puu) rms p rm s 

-0.0025 

0.003 

noise 

noise 

-0.0041 

0.0014 

-0.0019 

0.0007 
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Microphone 



Figure 1 . — Schematic of coordinate 
system and microphone locations. 



Figure 2. — Microphone auto-correlation in Mj = 0.98 jet (a) measured 
by a far field microphone at 9p = 160°; (b) calculated by the 
application of the causality principle to the shear noise source <u 7 ; p ; > 
measured from a large number of points in the plume. Plot from 
Schaffar 9 (1979). 



Figure 3. — Photographs of the jet facilities; (a) 1” unheated jets, (b) 2” heated jets. 


NAS A/CR— 2005-2 13817 


12 




Figure 4. — Normalized cross-correlation between sound pressure fluctuations p 7 and (a), (d) air density p; (b), (e) 
puu; (c), (f) pvv fluctuations measured in M a = 0.91 and 1.4 jets. For (a), (b), and (c) laser probe was at x/D = 10 
and centerline; for (d) and (e) probe is at x/D = 12 and centerline; for (f) it is x/D = 6 and r/D = 0.48. The 
microphone was kept at 50D and at indicated polar angles. The 90° microphone data are shifted by 0.05. 
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Figure 5. — Mach number (acoustic) and polar angle 
dependence of peak correlation coefficients. The laser 
probe was moved at different axial positions along (al), 
(bl), (cl) lip line r/D = 0.5 and (a2), (b2), (c2) centerline 
r/D=0 of 3 different plumes of M a = 0.6, 0.9, 1.48 and 
respectively, T r = 2.27, 2.7, 2.27 jets. The different color 
plots are for different microphone polar angles as 
indicated in the color code. All microphones were kept 
on a 100D arc. 


Figure 6. — Temperature dependence of the peak <p; 
p 7 > correlation. M a = 0.9 and static temperature 
ratios T r are as indicated. For the left column plots 
laser probe was moved along lip line and for the 
right column, it was moved along the centerline. 
The different color plots are for different 
microphone polar angles as indicated in color code. 
All microphones were kept on a 100D arc. 
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(a) \! a = 1.4 (a)M a = 0.91 


Figure 7. — Distribution of normalized cross-spectrum (coherence) between turbulent density fluctuations and far 
field sound pressure fluctuations at indicated Strouhal frequencies for two different unheated jets. The 
microphone was fixed at r/D = 50, 0 = 150° while the laser probe was moved from point to point. 



Figure 8. — Comparison of coherent spectrum between flow density and far field sound pressure fluctuations from 3 
different nozzle configurations: (a) baseline circular, (b) 16 lobe rectangular and (c) 4-tab circular nozzles. Laser 
position was moved at indicated axial stations along centerline of unheated M a = 0.91 jet. Plots from different x 
positions are shifted by 0.025. (d) Noise spectra from the fixed microphone at x/ D eqv =50, 0 P = 150°. 
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Figure 9. — (a) Coherence spectra between flow density fluctuations and 
far field sound pressure fluctuations in tabbed jet. The probe volume 
was moved radialy at indicated stations at fixed x/D eqv = 3. The 
microphone was fixed at x/D=50, 9 P = 150°. (b) Schematic of the flow 
field from 4 tabbed nozzle. 
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Figure 10. — Normalized correlation between air density fluctuations and (a), (b) selected individual microphones; 
(c) beam-formed signal created from all microphones. The laser probe was placed at x/D = 8, r/D = 0 in M a = 
0.91 jet. (d) A “unwrapped” view of the phased array indicating microphone locations. Laser probe volume lies 
directly above the array center between microphones 23 and 29. The small arrows in parts (a) and (b) indicate 
free-space propagation time from probe volume to microphone. 
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Figure 1 1 . — Maximum value of correlation between sound fluctuations at 
fixed microphone (0 P = 150°, R/D = 50) and (a) puu and (b) p fluctuations 
from various radial positions at x/D = 10 inside unheated M a = 1.18 jet. 
Note that correlation values are not normalized. 




Figure 12. — Schematic of large and small scale turbulent structures and noise radiation. 
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